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components,  message  delay. 


12l  ATTIIACT  (Cmmrn  m  mmm  it  ■■■■■■»  m*  »»  mmt  — Nri 

This  paper  presents  a  model  to  calculate  the  reliability  of  communication 
networks  with  multimode  components.  Previous  research  on  network  reliability 
has  focused  on  models  in  which  each  component  may  be  in  one  of  two  mode  s, 
namely,  operative  or  failed.  In  reality,  a  component  may  undergo  degradation 
in  performance  before  a  complete  outage,  and  will  therefore  operate  in  more 
than  two  modes.  Traditional  network  reliability  measures,  such  as  the  probabi 
litythat  a  pair  of  nodes  is  connected,  are  not  meaningful  in  a  multimode  mode 
Therefore,  the  mean  message  delay  of  the  network  is  defined  as  the  performanc 
measure.  An  exact  calculation  of  this  reliability  measure  is  not  feasible  due 
to  the  large  number  of  network  states,  corresponding  to  network  components 
being  in  different  modes.  We  have  developed  an  approximation  method  to  calcu 
late  this  reliability  measure.  This  method  requires  us  to  work  with  the  stat 
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Akatraet 

This  paper  protanta  a  modal  to  ealculata  tha 
roHabUlty  of  communication  natworka  with  multimoda 
eomponanta.  Pravloua  raaaarch  on  natwork  raHabilKy  fiaa 
focuaad  on  modaia  in  which  aach  componam  may  ba  in 
ona  of  two  modaa,  namaly,  oparativa  or  failed.  In  reality,  a 
component  may  undergo  dagradatlona  in  parformanca 
before  a  compiata  outage,  and  wHI  tharafora  operate  in 
mora  than  two  modaa.  Traditional  natwork  raiiability 
maaauraa.  auch  aa  tha  probability  that  a  pair  of  nodaa  la 
connactad.  are  not  meaningful  in  a  multimoda  modal. 
TiMrafora.  tha  mean  maaaaga  delay  of  tha  natwork  ia 
defined  aa  tha  parformanca  maaaura.  An  exact 
calculation  of  thia  rallabiitty  maaaura  la  not  faaaibia  due 
to  tha  large  number  of  natwork  atataa.  corraaponding  to 
natwork  eomponanta  being  In  dHfarant  modaa.  We  have 
developed  an  approximation  method  to  calculate  thia 
reliability  maaaura.  Thia  method  raquiraa  ua  to  work  with 
tha  atataa  of  tha  natwork  In  order  of  dacreaaing 
probablUty.  An  algortthm  OflOER-M  ia  davalopad  to 
ganarata  thaaa  atataa  In  tha  proper  order. 


1;  Iwkraductlan 

In  analyzing  tha  parformanca  of  a  communication 
natwork  with  unrallabla  eomponanta.  it  la  uaually  aaaumad 
that  aach  fallura-prona  component  can  ba  in  ona  of  two 
modaa.  either  operative  or  failed.  In  tha  oparativa  mode, 
tha  component  can  handle  load  at  Ka  capacity,  while  in 

tha  failed  mode.  N  ia  not  avallabla  at  all  Natwork 
reliability  maaauraa  auch  aa  tha  probability  that  a  given 
pair  of  nodaa  are  connactad.  tha  probability  that  all  paira 
of  nodaa  are  connected,  tha  axpactad  number  of 
communicating  node  paira,  and  ao  on  may  than  ba 
dafmad.  Numaroua  aigonthma  have  been  propoaad  (aaa. 
a.g..  (11)  for  computing  tha  aforemamlonad  natwork 
raiiability  maaauraa.  However,  tha  computation  time  la 
anormoua  axcapt  for  amaM-aiza  natworka.  In  tact,  it  haa 


Hhia  raaaarch  ia  aupportad  in  part  by  the  Air  Force 
Office  of  Scientific  Raaaarch  under  Contract  No 
AK)SR-B4-0269 


bean  proved  that  tha  exact  computatlona  of  thaaa 
raiiability  maaauraa  are  NP-hard  (21  (101  Other 
roaaarchart  (3]  have  propoaad  an  approach  that 
anumarataa  a«  poaalMa  natwork  atataa  and  caicuiaiaa  tha 
quantitlaa  of  imaraat  in  aach  atata  The  raliabilitv  maaaura 
ot  mtaraat  la  than  calcuiatad  aa  a  weighted  (  by  tha 
probability  of  that  atata  )  average  of  thaaa  quantitlaa. 
Again,  thia  approach  ia  Impractical  in  ganaral  bacauaa  tha 
number  of  natwork  atataa  growa  axponantially.  l.a..  If 
there  are  n  unrallabla  eomponanta.  thara  are  2'*  atataa  for 
tha  network.  Wa  tharafora  aaak  an  approximation  method 
to  minimize  tha  computation  tima.  In  (01  an  approximation 
approach  la  propoaad.  To  find  a  reliability  maaaura.  only 
tha  moat  probabla  atataa  ot  tha  network  are  analyzed. 
Tha  idea  la  that  whan  tha  atataa  conaidarad  accoum  for  a 
large  fraction  of  tha  atata  apace  (in  tarma  of  probability), 
wo  can  got  a  good  approximation  to  thia  raiiability 
maaaura.  Thia  approach  ia  practical  bacauaa  aoma  of  tha 
natwork  atataa  may  have  vary  amall  probabllitlaa,  and  tha 
network  rarely  oparataa  at  thaaa  atataa. 

While  two-mode  component  modaia  have  racaivad 
a  groat  deal  of  attention,  in  roallty,  a  component  may 
oparata  in  ona  of  N  (N>3)  modaa.  l.a..  a  component  may 
undergo  dagradatlona  baforo  a  compiata  outage.  For 
oxampla.  If  a  radio  channel  la  being  lammod,  ona  may  try 
to  combat  thia  lamming  by  ualng  a  more  powerful  arror- 
corractlwg  coda,  with  a  corraaponding  reduction  In  tha 
affoctiva  data  rata  or  channel  capacity  (111  Tharafora.  the 
channel  arm  operate  In  multipia  modaa.  with  aach  moda 
corraaponding  to  a  diffaram  channel  capacity.  In  (51  El- 
Nawaihi.  at.  al.  conaMar  an  N-atata  ayatem  which  ia  a 
function  ot  N-atata  (-mode)  eomponanta  and  thay  ara 
Intaraatad  in  tha  atructura  of  tha  ayatem  function.  Thay 
hava  not  however.  Indicatad  how  ona  may  dahna  tha  N 
atataa  ot  tha  natwork  in  a  meaningful  way  in  tarma  of  tha 
N  modaa  of  each  link.  In  thia  paper,  a  multimoda  modal 
for  network  reliability  will  ba  preaanted  and  uaad  to 
compute  aoma  properly  choaen  reliability  maaauraa 


2.  A  MuWmada  Modal 

Wa  aaauma  a  point-to-point  iMtwork  in  which  only 
links  may  fail  or  dagrada  and  nodaa  ara  failaafa  The 
modaa  of  a  link  corraapond  to  different  communication 
capacitiaa  ot  that  link.  i.a..  aach  mode  corraaponda  to  a 
diffarant  capacity.  Tharafora.  whan  a  link  la  not  operating 
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■t  H»  lUH  CaMCity.  M  WMV  optfWU  M  •  l0«Mr.  d«0ratfM 

MTvte*  rat*  facn  motfa  of  tit*  componom  l«  aMociatao 
witlt  a  known  probaMlltv  Thar*  ar*  M  Hnkt  bi  tha  natwork 
and  N  mod**  for  aach  Hnk.  caHad  modaa  0.tX.  -N-l 
(Wa  can  handia  componanta  with  loaa  than  M  modai  by 
dafining  tha  probabtlltlaa  of  axtranaous  mod**  a*  taro) 
Tha  |th  moda  of  Hnk  I  la  danotad  c/  and  tha  aaaociatad 
probability  l*  danotad  P|^*  Pr(C|l).  W*  hirthor  aatuma  that 
all  Hnk  fallura*  ar*  Indapandont  Not*  that  thar*  i*  a 
conttraint  for  aach  i: 

N-l 

I  I-U . M 

(•0 

Olvan  tha  natwork  topology  and  tha  tat  of  probabllltiat 
lor  link  modas.  wa  want  to  dafin*  a  quantity  to  moatura 
ih*  raiiabillty  of  tha  natwork.  If  wo  want  tha  probability 
that  a  givan  pair  of  nod**  ar*  connactad  by  a  path  of 
cartain  capacity,  w*  may  trantform  thi*  probtam  to  a 
typical  probiom  In  tha  two-mod*  modal  by  calling  tha  link 
modas  which  hava  at  laatt  that  capacity  at  oparativ*  and 
othar  modat  a*  failad  Tharafor*  computing  tha 
probability  of  tha  axistanc*  of  such  a  path  may  b* 
obtained  by  tha  slgoiithmt  dascribad  in  111  Wa  tharafora 
teak  othar  parformanc*  maaturai  which  ar*  meaningful 
lor  multimod*  modals,  and  which  raquir*  solutions 
significantly  dtflarant  from  thot*  for  two-mod*  modals 
Thera  hava  bean  attampts  to  us*  the  probability 
distribution  functiem  of  tha  maximum  flow  batwaan  a 
givan  sourca-dastination  pair  as  tha  parformanc* 
maatur*  in  tha  multimod*  modal.  For  axampi*.  Evans 
(61  studied  this  maasur*  for  a  natwork  in  which  tha 
capacity  of  any  adg*  i*  a  finita-stat*  mtagar-valuad 
random  variable  Kulkarni  and  Adlakha  (81  contidarad  a 
aimilar  problem  on  (t.t)-planar  networks  with 
axpoflMtialiy  diatrlbutad  Hnk  capaeitias.  In  this  paper,  tha 
natwork  mean  masaag*  dalay  is  ehosan  as  tha  raiiabillty 
measure  tine*  It  I*  related  to  the  degradation  in  channel 
capacity  for  the  links  In  tha  network.  With  a  given  node- 
to-node  traffic  raquiramant.  w*  are  mtarested  in  the 
probability 'that  tha  natwork  can  support  the  traffic  with  a 
given  average  delay.  Whan  the  mean  meaaag*  delay  is 
greater  than  this  threshold  at  a  '•ivan  network  state,  we 
say  that  the  natwork  cann^  satisfy  the  traffic 
requirement*  at  that  state.  Hence  w*  need  to  compute 
the  probability  that  the  mean  network  massage  delay 
does  not  exceed  soma  flnit*  value  For  aach  state  of  the 
natwork  (charactaniad  by  a  sat  of  link  modas),  w*  can 
compute  the  mean  messag*  delay  T  using  Kleinrocfc’s 
modal  (71: 


T--  I  - 

where  y  is  the  sum  of  the  arrival  rates  between  all  pairs 
of  nodes  1/u  IS  the  mean  packet  length,  Xj  is  the  trafhc 
load  on  link  i  and  C,  is  tha  capacity  of  link  i  at  that  state 
W*  assume  that  tha  components  change  states  relatively 
slowly  so  that  tha  queueing  phenomena  at  these  hnk* 
reach  equilibrium  and  the  avaraga  messag*  delay  **t  a 


given  state  of  the  network'  la  meaningful 

Having  dafinad  natwork  messag*  delay  a*  the 
raliablHty  maasur*,  w*  need  an  afficlant  way  to  compute 
It.  Since  there  ar*  network  states  In  our  modal,  we 
decide  to  ua*  an  approach  yyhich  enumaratos  only  the 
most  probaM*  states  and  analyiat  the  parformanc*  for 
that*  state*.  W*  hava  davalopad  an  algorithm  to  gonarat* 
natwork  atatas  in  order  of  docraasing  probability  for  our 
multimod*  modal. 


8.  AlgaiHItm  ONOfR-M 

There  ara  M  Hnk*  and  N  modat  for  aach  link  In  the 
natwork  The  modat  are  renamed  such  that  the  mod* 
probabilities  ar*  In  decreasing  order,  i.*.,  P|*^>Pi^>'  ,  ■> 
Pj^*^  for  all  I.  Each  natwork  state  S  is  rapresentad  by  an 

M-v*ctor.  S*()y,i2 . ImI-  whar*  |j  denotes  the  mod*  of 

Hnk  I.  Tha  most  probable  state  is  (0,0..  .  .,0)  with 
probability  W*  want  to  gonarat*  the  m  most 

probable  states  in  order  of  decreasing  probabilities  Our 
approach  I*  *  modification  of  the  algorithm  ORDER  (9l 
which  It  dasigned  for  two-mod*  natwork*. 

Lot  Q>MMN-1).  For  pj.  15  i<  M.  1<  (<  N-1.  w* 
arrange  that*  Q  parameters  in  dacraasing  order  and  us* 
01-  -Qq  denote  this  ordered  saquonca  so  that  each 
Ol^  corraaponds  to  a  diffarant  p/  Not*  that  p^°.  ..p^ 

ar*  not  raquirad  in  tha  algorithm  and  ar*  not  included  in 
this  set  of  Q  paramatars  Wa  first  dascrib*  the  way  new 
states  ar*  generated.  Let  A*{S,.S|^  .  .  .S^)  be  an  ordered 
sat  of  natwork  ttatas  such  that  Pr(S,)>.  .  .>Pr(S^).  Define 
an  operation  raplac*  (i.i).  danotad  r]^.  as  follows:  R|<Sf 
mean*  that  tlta  ith  aiamant  of  state  S,  Is  changed  to 
mod*  i  and  the  othar  aiamantt  remain  tha  same  The 
probabilify  of  Ih*  new  atat*  ganaratad  is  easy  to 
compute: 

PrfRJSfHPfiS,!*  pJ/PrlS^i)l  (II 

where  S^l)  I*  the  ith  element  of  This  raplac*  operation 
may  also  be  defined  on  the  whole  sat  of  states  Thus 
B>R|iA>{R|is,,.  .  .Rj'S^}  Not*  that  the  size  of  set  B  may 
be  smaller  than  that  of  sot  A  because  two  states  S,  and 
of  A  may  differ  only  In  elemant  i  and  so  ar*  changed 
to  the  same  state  by  R,J  In  order  to  avoid  ganarating 
duplicate  new  states  by  thi*  raplac*  operation,  wa  need 
to  record  the  identities  of  those  comportants  whose 
modas  have  already  been  changed  in  previous  steps  W* 
us*  a  sot  l|^  to  record  the  components  used  to  genarate 
new  states  thus  tar  ff  the  ith  component  has  been  used 
before,  then  only  thot*  states  in  A  whose  ith  element  is  0 
Will  be  used  to  genarate  new  states  in  RjlA.  Define  an 
operation  insert,  denoted  B«  A.  as  follows:  C'B"  A  •the 
ordered  sat  which  results  when  every  state  of  B  is 
inserted  into  the  ordered  tat  A.  Define  an  operation 
salect.  denoted  T(A).  as  follows:  T(A)*tha  ordarad  set 
which  contains  tha  first  m  states  of  the  ordarad  s*< 
A  The  algorithm  consists  of  two  phases  and  Q  stages 
with  stag*  k  corresponding  to  parameter  Wa  star 
with  state  S^*(0,0..  .0)  To  ganarst*  new  states  we 


3 


MCCMsivcty  take  4^  q2.  .Vq  <***  corratpondtng 
lt}  %  and  raplaca  tha  ttk  aiatnam  of  pravtoutiy  ganaratad 
atataa  by  Wttan  tfia  total  numbar  of  atatat  ganaratad  la 
at  laatt  m.  wa  go  to  pbaaa  2.  in  pliaaa  2.  wa  ganarata 
now  atataa  by  a  abnliar  procadura.  but  wUi  ratain  only  tha 
m  moat  probaWa  atataa. 

Wa  now  glva  tha  aigonthm  OftOEh-M: 

Phaaa  i: 

I  mitlallta:  S,-<0.  ,0).  Ao*{S^}.  lo*4 

2.  For  k-14 . m-l.  rapaat  3.4.5,6 

3.  Find  tha  P)^  corraapondmg  to  4|^ 

If  iclifi.  than  flag«l.  alaa  flag«0 

4.  For  r-l. .  .  .4Aj.||.  do: 

it  flag'1  and  S^i)4  0.  than 

no  naw  atata  ia  ganaratad  from  S^: 

alaa  ganarata  naw  atata  ftj'Sr  and  computa 
ita  probability 

(Tha  aat  of  naw  atataa  ganaratad  in  thia  stap 
■a  danotad 

5.  A|^.|  and  ia“ia-t^^*^' 

6.  If  tA|J>m.  go  to  phaaa  2 
Phaaa  2: 

y.  Initialita.  A*T(A|^)  whara  L  is  tha  loop  indax  k 
whan  wa  laava  phaaa  1  Tha  naw  aat  A  is 
than  rartamad  A^^. 

8.  For  k>L*1 . Ol  rapaat  9.10.11.12 

9  Taka  and  tha  corraaponding  Pjl 
If  <ci|[.^.  than  flag'1.  alaa  fltp'O. 

10  For  r-l.  do. 

If  flag-t  and  S^i)n  0.  than 

no  naw  atata  is  ganaratad  from  S,.. 

alsa  ganarata  naw  atata  and  computa 
its  probability 

(This  stop  may  ba  danotad  B|^.f*R,lA^.f ) 

II  A|^_  ^  and  Ifc^ly-t^Ci) 

12  Ak-T(A) 


Aq  contains,  m  dacraasing  oroar.  tha  m  most 

probabis  ttatas. 

Tha  following  obaarvation  on  ORDER-M  can  ba 

mada. 

Obaarvation  1;  Tha  oparations  Rj^.  1<i<M.  l<i<N-1. 
ara  axhaustiva  of  tha  ways  naw  sutas  can  ba  ganaratad 
from  old  onas. 

Lamma  1:  In  aach  staga.  for  any  two  statas  S,  and 
S„  which  hava  alraady  baan  ganaratad  such  that  Pr(S,]> 
PrtS„l  wa  must  hava  Pr(R)*S,l>Pr(Ri*S„]  H  both  of  thasa 
naw  statas  axist 

Proof:  Whan  both  RjlS^  and  Rj^Sp  axist  wa  mutt 
hava  S^i)*Sp(l)>0.  It  follows  from  Equation  (1)  that 
PfiRjlS^>PrlR|ISnl 

Lamma  1  statas  that  tha  ttatas  in  B|^_^  ganaratad  in 
staps  4  and  10  ara  in  ordar  of  dacraasing  probability.  It 
alto  aiturat  that  statas  laft  out  by  tha  salact  oparation 
T(A)  will  not  glva  rita  to  now  ttatas  in  tha  top  m  spots  m 
tha  naxt  staga. 

Tha  naxt  thaoram  it  a  diract  consoquanco  of  tha 
abova  obaarvation  and  lamma. 

Thaoram  I:  At  tha  and  of  ORDER-M,  tha  m  most 
probabla  statas  ara  ganaratad. 

Wa  now  astimata  tha  tima  complaxiry  of  this 
algorithm  For  aach  k  in  phase  2.  Rj^A  takas  tima  0(m* 
M)  and  chocking  whathar  icl|t  Ukat  tima  0(M).  In  tha 
Intartion  wa  can  uta  a  sorting  algorithm  with  tima 
complaxity  0(m*  logm).  Tha  salact  oparation  can  bo 
porformad  during  tha  Insartlon.  i.o..  In  tha  sama  toning 
routine.  The  complaxHy  analysis  of  a  staga  in  Phasa  1  is 
similar  to  Phasa  2.  However,  instead  of  working  with  a  list 
of  m  ttatas  (  the  m  most  probable  states  ganaratad  thus 
far  ).  one  it  working  with  m'<m  states  in  phase  1.  Thera 
ara  a  total  of  Q  raplaca  operators  in  both  phases.  So  tha 
total  time  in  tha  worst  case  is  0(0  m  M'rQ  m  logm) 

4.  Covorage  of  ttia  State  Space 

In  order  to  astimata  how  many  statas  must  bo 
ganaratad  to  achiavo  a  cenain  covaraga  of  tha  state 
space,  wa  consider  tha  case  In  which  for  any 

mode  j  and  links  i  and  k.  To  simplify  tha  notation,  wo  let 
pt>P|l  for  all  i.j  and  assume  as  before  that  tha  labeling  tor 
tha  statas  is  such  that  P^>P^>.  .>p^*V  Each  state  of 

the  network  It  reprasantad  by  a  random  vector 

whara  Xj  is  tha  mode  of  the  ith  link.  0<X^N-1.  1< 
i<M  Note  that  Xj's  ara  indopandent  random  variablas  in 
our  modal  Define  a  random  variablo  X*  J,.,X,.  Tha  most 

probabla  state  ia  (0.0 . 0|  with  X>0.  tha  naxt  most 

probabla  has  X*1.  tha  third  most  probable  may  hava  X*1 
or  2.  and  to  on.  Tha  most  probable  system  ttatas 
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by  our  tlgorttftin  cotrotpond  to  tho  oyont  (X< 
A),  whoro  A  It  •  pottttvo  intogor  Th«  Cltornoff  bound 
I4.  121  can  bo  utod  to  dorivo  an  upoar  bound  on  tna 
probabWty  of  tba  tat  of  uncovorad  ttatat.  ta..  PrfX>Ai. 
Wo  firtt  introduca  toma  notatlont  in  tfta  following 
paragraph. 

Tha  momant-Mnorating  function  of  X|.  i~t . M.  it 

♦(t)-EIa>ip<iX,)l-jJ^  p*ai«p(ti)  Tha  tami-invartant 
ganarating  function  it  p(t>*lnp(s).  For  random  vanabla  K 
tha  momant-ganaratlng  function  it 
Itt  aami-invariant  ganarating  function  it  ^X^*l  * 

M  ln4(t).  Tha  Charnoff  bound  ttatat  tha  following 

Pr(X>A)<axp(Px<i)-iA)-aiip(Mw(t)-iAl  (21 

whtra  tha  optimum  valuo  of  t  (which  givai  tha  tightoit 
bound  of  tha  axponontial  form)  it  talactad  in  accordanco 
with 


•.  Oalay  Anaiyalt  of  a  Natarark  wfoi  MuHlatada  Unka 
Wo  now  apply  OftOCR-M  to  ttM  raHabNfty  anatytit  of 
tha  natwork  In  Fig.  1.  R  hat  S  nodat.  MM  Unka  and  aach 
link  hat  M>5  modat  to  0*24  Tha  axtomal  traffic 
raquiramontt  batwaan  nodo  paira  aro  atiumad  to  ba 
Poltton  with  maan  Td-  whara  t  It  tha  tourca  and  t  It 
tha  dattinatlon.  t.t«i’.2..  .  .5.  Tha  ttatat  and  tha 
corrotponding  capacitiat  of  tha  links  art  givon  in  Tsbia  V 
Tha  attociatad  probsbilitiot  art  ghran  in  Tabla  2 

i 


0 

1 

2 

3 

4 

1 

50 

40 

30 

10 

20 

2 

50 

40 

30 

10 

20 

3 

50 

40 

30 

10 

20 

4 

50 

30 

40 

20 

10 

5 

50 

30 

40 

20 

10 

6 

SO 

30 

40 

20 

10 

dtiltl 

M - 


«• 


•A 


Difforontiating  tha  laft-hand  tida  glvtt 


(3) 


1 


or 


(41 


So  talacting  tha  optimum  t  involvas  toiving  a  polynomial 
aquation  of  dagraa  N-1  in  a* 


In  applying  this  bound  to  a  ganaral  N-moda 
proWom.  thara  art  soma  complicationt.  howavar.  Whan 
N>3.  tha  cata  of  two  componantt  at  moda  1  contributat 
tlia  tama  valba  to  X  as  tha  cata  of  ona  componant  at 
moda  2.  Thata  two  catas  may  hava  vary  diffarant 
probabilltttias.  Such  complications  will  not  arisa  in  a 
two-moda  modal,  whara  tha  Charnoff  bound  takas  tha 
following  form  Pr<X>  A\<  axplM  H(A/M.  1  *A/I4HM-A1 
lnp®*A  Inp’l  whara  H(q.1-q}— qlnq-(1-qlln(l-ql  and 
p^tp^al.  In  this  cata.  for  a  givan  covarsga  of  tha  stata 
spaca.  say  95%.  wa  can  saiact  A  by  using  tha  abova 
formula  for  Pr(X>iy^0.05  Tha  numbar  of  ttatat  noadad 
Is  than  m- This  givas  a  quick  way  o» 
datarmining  tha  vaiiia  m  for  our  algorithm  Howavar.  sinca 
tha  Charnoff  bound  is  an  inaquality  tha  rasult  is 
pastimittic.  masning  that  m  it  usually  largar  than  sctuslly 
naadad.  although  it  is  asymptotically  tha  bast  among  all 
axponantial  forms  (121 


Tabla  1.  Cji  in  kbps  15156.  05)54 
I 


0 

1 

2 

3 

4 

1 

99 

006 

002 

.001 

001 

2 

95 

.025 

01 

.01 

005 

3 

99 

.007 

001 

.001 

.001 

4 

97 

018 

0054 

.0036 

.003 

5 

95 

.03 

01 

.005 

005 

6 

97 

024 

003 

.0018 

0012 

Tabla  2. 

Pi*  . 

151  <6.  05154 

Wa  first  arranga  tha  24'probabtlltiat  Pj’,  15156.  1515 
4  in  dacraaslng  ordar  and  rocord  tha  (i.i)  pair  of  tha 
corrotponding  q^  by  a  24x2  matrix 

i  I 

1  5  1 

2  2  1 

3  6  1 

4  4  1 

5  2  2 

6  5  2 

7  2  3 

8  3  1 

9  1  1 

10  4  2 

11  5  3 

k  12  2  4 

13  5  4 

14  4  3 

15  6  2 

16  4  4 

17  1  2 

18  6  3 

19  6  4 

20  3  2 

21  1  3 


22 

23 

24 
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3  3 

1  4 

3  4 

Table  3  b),.  k*l.2 . 24 

t.1.  Mm  Estimata 

To  analyxa  tita  mean  masMga  delay  at  a  raliabilltv 
maatura  of  tba  network,  algorithm  OftOER-M  it  uted  to 
generate  the  m  mott  probable  itataa.  For  each  ttata  of 
tha  network,  we  calculate  the  network  mettage  delay  by 
Klainrock't  modal  uting  a  fined  routing  table.  Define  an 
Indicator  variable  f^  for  each  itate  S|  at  follows: 

f|a  1  If  the  network  delay  it  <200  mt  in  Sj  (6) 

0  otherwise 

then 

F*Pr(network  mean  mettage  delay  it  <200  mt) 


■[PrfSi)fj 

Hera  tha  minimum  delay  among  all  states  it  106  7 
ms  (this  value  is  achieved  at  ttata  S^*(0.0.0.0.0.0))  and  tha 
value  200  mt  it  chosen  as  an  acceptable  value  for  the 
network  delay,  if  the  above  summation  covers  all  states, 
then  we  gat  tha  enact  value  of  F  It  it  alto  possible  to 
compute  lower  and  upper  bounds  for  F  Having  calculated 
f^  lor  the  mott  probable  states,  if  we  let  fj*0  for  all  othar 
states,  wo  get  a  lower  bound:  H  wo  let  fj*!  for  all  othar 
states,  we  get  an  upper  bound. 

•The  following  routing  metrin  it  uted  for  this 
enample 


destination 

1  2 

3 

4 

5 

1 

1.44 

2 

2.3.5 

3 

3  2.1 

3.5 

4 

4.5 

S 

5.3.2 

5.3 

5.4 

Tha  matrin  entries  represent  tha  vamcas  along  which  the 
data  packets  are  forwarded. 

We  generate  the  m«20  mott  probable  states  for  this 
network  by  ORDER-M  These  states  and  their  associated 
probebilitiat  tra  shown  in  Table  4  The  sum  of  the 
probabilities  of  these  20  states  is  calculated  to  be 
0  98396  The  reliability  measure  F.  tha  probability  that  the 
rtetwork  satisfies  tha  antamel  traffic  requirements  (  le. 
tha  network  delay  does  not  onceed  200  ms  I.  is 
tppronimated  as  0  92477  using  these  20  stales  The  lower 
bound  and  upper  bound  for  the  value  F  as  a  function  of 
m  the  number  of  states  considered,  are  plotlad  in  Fig  2 


for  i<m<20.  We  have  thus  iNuatrated  the  enaiysis  of  a 
network  with  multimode  components  and  presented  an 
Indication  of  how  reliable  tha  network  la  In  terms  of 
providing  a  apecified  level  of  service 


rank 

atate  vector 

probebUlty 

1 

(0.0,0.0.0.0) 

832264 

2 

(0.0.0.0.1.0) 

.026282 

3 

(0.1.0.04.0) 

.021902 

4 

(0.0.0.0.0.1) 

.020502 

5 

(0.0,0.1.0.0) 

.015444 

6 

(0.2.0.0.0.0) 

.008761 

7 

(0.0.0.04.0) 

.008761 

8 

(0.3.0.0.0.0) 

.008761 

9 

(0.0.1. 0,0.0) 

.005885 

10 

(1.0,0.0.0.0) 

005044 

11 

(0.0.04,0.0) 

.004633 

12 

(0.0.0.04.0) 

.004380 

13 

(0.4.0.0.0,0) 

.004380 

14 

(0.0.0.0,4.0) 

.004380 

15 

(0.0.04.0.0) 

.003089 

16 

(0.0,0.00.2) 

002574 

17 

(0.0.0.4.0.0) 

002574 

18 

(2.0.0.0.0.0) 

001681 

19 

(0.0.0.0.0.3) 

001544 

20 

(0.0.0.0.0.4) 

001030 

Table  4  The  20  mott  probable  states 
and  their  probabilities 


$J.  Eadmatlns  the  Oiatributlen  Function 

An  alternate  approach  to  tha  delay  analysis  for  our 
model  is  to  find  the  probabllrty  distribution  function  of  the 
network  mean  message  delay.  Define  F(t>>Frtnatwork 
mean  message  delay  is  <i  ms)  and  an  indicator  variable 
f,(t)  for  state  S|: 

tj(t)>1  If  tha  network  delay  It  <t  ms  In  S|  (7) 

0  otherwlsa 

thenF(t)  -[,l*rtS,)fj(t) 

The  exact  determination  of  F(t)  requires  tha 
enumeration  of  the  whole  state  space  However,  whan  we 
generate  tha  m  mott  probable  states,  we  can  calculate 
F„,|t)>^”i  PTtSi)fi(t).  This  value  Is  than  a  lower  bound  on 
F|t).  I.e..  F(t)>F^(t)  for  all  t.  For  our  example.  F^(t)  using 
m«20  and  tha  mott  probable  states  in  Tabia  4  it 
calculated  and  plotted  in  Fig  3 


E.  Conclusions 

In  this  paper,  we  have  prasantad  a  multimode  model 
for  calculating  network  reliability  and  define  the 
probability  that  the  traffic  requirements  are  satisfied  as 
the  reliability  measure  An  algorithm  OROER-M  it 
developed  which  enumerates  network  states  m  order  of 
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(McrMstng  probability  W«  tlion  apply  It  to  cotnputa  tha 
raliabitity  dobnoo  abova. 

Uiinp  our  alfloiitMn.  wa  ara  abla  to  eompula  tha 
probability  diatribution  function  of  tha  natwork  maan 
masaapa  daiay  to  a  canain  accuracy  m  a  vary  afficiant 
mannar  Tha  advantaga  of  our  modal  is  that  tha  iraHic 
raqutramant  of  tha  natwork  la  takan  into  account  in  tha 
raliabiliiy  analysis 


V" 

mmms  ana  nbn*  dnaKnhM 

SI  WOSTSMB 

lb  •  mm  moat  irnmimi  ■rn/Mcnr 

B>s  ikns 

figura  1  Exampla  Communication  Natwork 


Figura  2:  Uppar  and  Lowar  Sounds  tor  f 
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Figura  3  An  Estimata  ol  tha  Probability  Distribution 
Function 


